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Recently, there has been growing interest in using adiabatic quantum computation as an architec- 
ture for experimentally realizable quantum computers. One of the reasons for this is the idea that 
the energy gap should provide some inherent resistance to noise. It is now known that universal 
quantum computation can be achieved adiabatically using 2-local Hamiltonians. The energy gap in 
these Hamiltonians scales as an inverse polynomial in the problem size. Here we present stabilizer 
codes which can be used to produce a constant energy gap against 1-local and 2-local noise. The 
corresponding fault-tolerant universal Hamiltonians are 4-local and 6-local respectively, which is the 
optimal result achievable within this framework. 

PACS numbers: 03.67.Pp 



Adiabatic quantum computation was originally pro- 
posed by Farhi et al. as a method for solving combi- 
natorial optimization problems pj. In this scheme, one 
starts with a Hamiltonian whose ground state is easy to 
construct, and gradually varies the Hamiltonian into one 
whose ground state encodes the solution to a computa- 
tional problem. By the adiabatic theorem, the system 
will remain in the instantaneous ground state provided 
that the Hamiltonian is varied sufficiently slowly. More 
precisely, any closed system acted on by H(t/T) from 
t = to T will remain in the ground state with high 
probability provided that T is sufficiently large. Diffcr- 
ent formulations0, of the adiabatic theorem yield 

different conditions on T, but essentially the minimal T 
scales polynomially with the inverse eigenvalue gap be- 
tween the ground state and first excited state. 

Recently, there has been growing interest in using adi- 
abatic quantum computation as an architecture for ex- 
perimentally realizable quantum computers. Aharonov 
et a/.[||, building on ideas by FeynmanQ and KitaevQ, 
showed that any quantum circuit can be simulated by 
an adiabatic quantum algorithm. The energy gap for 
this algorithm scales as an inverse polynomial in G, the 
number of gates in the original quantum circuit. G is 
identified as the running time of the original circuit. By 
the adiabatic theorem, the running time of the adiabatic 
simulation is polynomial in G. Because the slowdown 
is only polynomial, adiabatic quantum computation is a 
form of universal quantum computation. 

Most experimentally realizable Hamiltonians involve 
only few-body interactions. Thus theoretical models of 
quantum computation are usually restricted to involve 
interactions between at most some constant number of 
qubits k. Any Hamiltonian on n qubits can be expressed 
as a linear combination of terms, each of which is a tensor 
product of n Pauli matrices, where we include the 2x2 
identity as a fourth Pauli matrix. If each of these tensor 
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produets contains at most k Pauli matrices not equal to 
the identity then the Hamiltonian is said to be fc-local. 
The Hamiltonian used in the universality construction of 
[f| is 3-local throughout the time evolution. Kempe et 
al. subsequently improved this to 2-local in Q. 

Schròdinger's equation shows that, for any constant g, 
gH(gt) yields the same time evolution from time to 
T/g that H(t) yields from to T. Thus, the running 
time of an adiabatic algorithm would not appear to be 
wcll defined. However, in any experimental realization 
there will be a limit to the magnitude of the fields and 
couplings. Thus it is reasonable to limit the norm of each 
term in H(t). Such a restriction enables one to make 
statements about how the running time of an adiabatic 
algorithm scales with some measure of the problem size, 
such as G. 

One of the reasons for interest in adiabatic quantum 
computation as an architecture is the idea that adiabatic 
quantum c omp uters may have some inherent fault tol- 
erance 0, 0, 0, Q] . Because the final state de- 
pends only on the final Hamiltonian, adiabatic quantum 
computation may be resistant to slowly varying control 
errors, which cause H{t) to vary from its intended path, 
as long as the final Hamiltonian is correct. An exception 
to this would oceur if the modified path has an energy 
gap small enough to violate the adiabatic condition. Un- 
fortunately, it is generally quite difficult to evaluate the 
energy gap of arbitrary local Hamiltonians. 

Another reason to expect that adiabatic quantum com- 
putations may be inherently fault tolerant is that the 
energy gap should provide some inherent resistance to 
noise caused by stray couplings to the environment. In- 
tuitively, the system will be unlikely to get excited out of 
its ground state if kbT is less than the energy gap. Un- 
fortunately, in most proposed applications of adiabatic 
quantum computation, the energy gap scales as an in- 
verse polynomial in the problem size. Such a gap only 
affords protection if the temperature scales the same way. 
However, a temperature which shrinks polynomially with 
the problem size may be hard to achieve experimentally. 

To address this problem, we propose taking advantage 
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of the possibility that the decoherence will act indepen- 
dently on the qubits. The rate of decoherence should 
thus depend on the energy gap against local noise. We 
construct a class of stabilizer codes such that encoded 
Hamiltonians are guaranteed to have a constant energy 
gap against single-qubit excitations. These stabilizer 
codes are designed so that adiabatic quantum compu- 
tation with 4-local Hamiltonians is universal for the en- 
coded states. We illustrate the usefulness of these codes 
for reducing decoherence using a noise model, proposed 
in [ïoj . in which each qubit independently couples to a 
photon bath. 

To protect against decoherence we wish to create an 
energy gap against single-qubit disturbanecs. To do this 
we use a quantum error correcting code such that apply- 
ing a single Pauli operator to any qubit in a codeword 
will send this state outside of the codespace. Then we 
add an extra term to the Hamiltonian which gives an en- 
ergy penalty to all states outside the codespace. Since 
we are only interested in creating an energy penalty for 
states outside the codespace, only the fact that an error 
has oceurred needs to be detectable. Since we are not 
actively correcting errors, it is not necessary for distinct 
errors to be distinguishable. In this sense, our code is 
not truly an error correcting code but rather an error de- 
tecting code. Such passive error correction is similar in 
spirit to ideas suggested for the circuit model in [ÏÜj . 

It is straightforward to verify that the 4-qubit code 

\0l) = -(|0000) + i|0011) + i|1100) + |llll)) (1) 
\1 L ) = ^ (-10101) +t|0110> +*|1001) — |1010» (2) 
satisfies the error-detection requirements, namely 

(OlMOl) = = (0 L \a\l L ) = (3) 

where o is any of the three Pauli operators acting on one 
qubit. Furthcrmore, the following 2-local operations act 
as encoded Pauli X, Y, and Z operators. 

X L = Y®I®Y®I 

Y L = -I®X®X<g>I (4) 

Z L = Z®Z®I®I 

That is, 

X l \0l) = X l \1 l ) = \0 L ), 

Y L \0 L ) = i\l L ), Y L \1 L ) = -i\0 L ), 
Z L \0 L ) = \0l), Z L \l L ) = -\\ L ). 

An arbitrary state of a single qubit a\0) +/3|1) is encoded 
as a\0 L )+P\l L ). 

Starting with an arbitrary 2-local Hamiltonian H on 
N bits, we obtain a new fault tolerant Hamiltonian on 
AN bits by the following procedure. An arbitrary 2-local 
Hamiltonian can be written as a sum of tensor prod- 
uets of pairs of Pauli matrices acting on different qubits. 



After writing out H in this way, make the following re- 
placements 

X->X L , Y^Y L , Z->Z L 

to obtain a new 4-local Hamiltonian Hsl acting on AN 
qubits. The total fault tolerant Hamiltonian Hs is 

H s = H SL + H SP (5) 

where Hsp is a sum of penalty terms, one acting on each 
encoded qubit, providing an energy penalty of at least E p 
for going outside the code space. We use the subscript S 
to indicate that the Hamiltonian acts on the system, as 
opposed to the environment, which we introduce later. 
Note that Hsl and Hsp commute, and thus they share 
a set of simultaneous eigenstates. 

If the ground space of H is spanned by \ip^) ■ ■ . \ip( m ') 
then the ground space of Hs is spanned by the encoded 

states IVx ) ■ ■ • I^Pl" 1 ^)- Furthermore, the penalty terms 
provide an energy gap against 1-local noise which does 
not shrink as the size of the computation grows. 

The code described by equations and [21 can be ob- 
tained using the stabilizer formalism [ÏQ, UJX- In this 
formalism, a quantum code is not described by explicitly 
specifying a set of basis states for the code space. Rather, 
one specifies the generators of the stabilizer group for 
the codespace. Let G n be the Pauli group on n qubits 
[i.e. the set of all tensor produets of n Pauli operators 
with coefficients of ±1 or ±i). The stabilizer group of a 
codespace C is the subgroup S of G n such that x\ip) = \tp) 
for any x € S and any \ip) € C. 

A 2 fc dimensional codespace over n bits can be specificd 
by choosing n — k independent commuting generators for 
the stabilizer group S. By independent we mean that no 
generator can be expressed as a product of others. In our 
case we are encoding a single qubit using 4 qubits, thus 
k = 1 and n = 4, and we need 3 independent commuting 
generators for S. 

To satisfy the orthogonality conditions, listed in equa- 
tion[31 which are necessary for error detection, it sufhces 
for each Pauli operator on a given qubit to anticommute 
with at least one of the generators of the stabilizer group. 
The generators 

gi = X <g> X ® X ® X 

g 2 = Z (g) Z ® Z ® Z 

g 3 = X®Y ®Z®I (6) 

satisfy these conditions, and generate the stabilizer group 
for the code given in equations and [3 
Adding one term of the form 

H p = -E p ( gi +g 2 + g 3 ) (7) 

to the encoded Hamiltonian for each encoded qubit yields 
an energy penalty of at least E p for any state outside the 
codespace. 

2-local encoded operations are optimal. None of the 
encoded operations can be made 1-local, because they 
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would then have the same form as the errors we are trying 
to detect and penalize. Such an operation would not 
commute with all of the generators. 

Intuitively, one expects that providing an energy gap 
against a Pauli operator applied to any qubit protects 
against 1-local noise. We illustrate this using a model of 
decoherence proposed in |ïcj . In this model, the quantum 
computer is a set of spin-1/2 partides weakly coupled to 
a large photon bath. The Hamiltonian for the combincd 
system is 



H = H S + H E 



XV, 



where Hs(t) is the adiabatic Hamiltonian that imple- 
ments the algorithm by acting only on the spins, He is 
the Hamiltonian which acts only on the photon bath, and 
XV is a weak coupling between the spins and the photon 
bath. Spccifically, V is assumed to take the form 



V = J2 



dus 



g(w)a u <r^ + g* {u)a1<j { ^ 



where a± are raising and lowering operators for the ith 
spin, a^, is the annihilation operator for the photon mode 
with frequency u>, and g(ui) is the spectral density. 

From this premise Childs et al. obtain the following 
master equation 



dp 



(8) 



a . b 



where 



M ab = ^2[N ba \g ba \ 2 (a\a^\b)(b\a^\a) 

i 

+ (A r a h + l)|5a6| 2 (&|^ ) |a)(a|4 ) l ò > 

is a scalar, 

£ ab (p) = \a)(a\p + p\a)(a\-2\b)(a\p\a)(b\ 



is an operator, \a) is the instantaneous eigenstate of H$ 
with energy u a , 



N, 



1 



ba 



exp[(3(uj b —Lü a )]—1 
is the Bose-Einstein distribution at temperature 1 //?, and 



9ba 



_ / Xg{uj b - u a ), 
0, 



w& > iü a , 

oj b < u> a . 



(9) 



Suppose that we encode the original iV-qubit Hamil- 
tonian as a 4iV-qubit Hamiltonian as described above. 
As stated in equation [5] the total spin Hamiltonian Hg 
on AN spins consists of the encoded version Hsl of the 
original Hamiltonian H$ plus the penalty terms Hsp. 



Most adiabatic quantum computations use an initial 
Hamiltonian with an eigenvalue gap of order unity, in- 
dependent of problem size. In such cases, a nearly pure 
initial state can be achieved at constant temperature. 
Therefore, we'll make the approximation that the spins 
start in the pure ground state of the initial Hamiltonian, 
which we'll denote |0). Then we can use equation[H]to ex- 
amine dp/dt at t = 0. Since the initial state is p = |0)(0|, 
£ ab {p) is zero unless \a) — |0). The master equation at 
t = is therefore 



dp 
dt 



t=o 



= -i[H s ,p}-J2M ob £ ob {p). (10) 



Hsp is given by a sum of terms of the formQ and it 
commutes with Hsl- Thus, Hs and Hsp share a com- 
plete set of simultaneous eigenstates. The eigenstates of 
Hs can thus be separated into those which are in the 
codespace C (í.e. the ground space of Hsp) and those 
which are in the orthogonal space C . The ground state 
0) is in the codespace. M ob will be zero unless |6) G C ± 7 
because a± = (X ± iY)/2, and any Pauli operator ap- 
plied to a single bit takes us from C to C . EauationllQI 
therefore becomes 



dp 
dt 



t=o 



i[H s ,p]+ M ob£ob{p) (11) 

hec 1 - 



Sincc |0) is the ground state, Lu b > u>q, thus equation 
EDshows that the terms in M ob proportional to |<7of>| 2 wm 
vanish, leaving only 

M ob =^iV i , |56o| 2 (0|ai i) |6><ò|4 i) |0). 

i 

Now let's examine A^bo- 

u b - uj = (b\(H SL + H S p)\b) - {0\(H SL + H S p)\0). 
0) is in the ground space of Hsl, thus 

(b\H SL \b) - <0|-ff Si |0) > 0, 

and so 

uj b -Lü Q > (b\H S p\b)-{0\H S p\0). 
Since \b) £ C 1 - and |0) £ C, 

(b\H SP \b) - (0\H SP \0) = E p , 

thuS Lú b — Lüq > Ep. 

A sufnciently large (3E p will make N ba small enough 
that the term X^ec- 1 - Mo b £(p) can be neglected from the 
master equation, leaving 



dp 
dt 



t=o 



which is just Schròdinger's equation with a Hamiltonian 
equal to Hs and no decoherence. Note that the preceding 



4 



derivation did not depend on the fact that a± are raising 
and lowering operators, but only on the fact that they 
act on a single qubit and can therefore be expressed as a 
linear combination of Pauli operators. 

Nbo is small but nonzero. Thus, after a sufficicntly 
long time, the matrix elements of p involving states other 
than |0) will become non-negligible and the preceding 
picture will break down. How long the computation can 
be run before this happens depends on the magnitude of 
X^bec- 1 - M b£(p), which shrinks exponentially with E p /T 
and grows only polynomially with the number of qubits 
N. Thus it should be sufhcient for 1/T to grow logarith- 
mically with the problem size. In contrast, one expects 
that if the Hamiltonian had only an inverse polynomial 
gap against 1-local noise, the temperature would need to 
shrink polynomially rather than logarithmically. 

Now that we know how to obtain a constant gap 
against 1-local noise, we may ask whether the same is 
possible for 2-local noise. To accomplish this we need to 
find a stabilizer group such that any pair of Pauli op- 
erators on two bits anticommutes with at least one of 
the generators. This is exactly the property satisfied by 
the Standard^?] 5-qubit stabilizer code, whose stabilizer 
group is generated by 

gi = X®Z®Z®X®I 

g 2 = I®X®Z®Z®X 

g 3 = X®I®X®Z®Z 

.94 = Z®X®I®X®Z. (12) 

The codewords for this code are 

|0 L ) = - [ 100000} + |10010) + |01001) + |10100) 

+ |01010) - |11011) - |00110) - |11000) 
— 111101) - |00011) - |11110) - |01111) 
— 110001) - |01100) - |10111) + |00101) ] 

\i L ) = j [ |imi) + loiioi) + 110110} + |oioii) 

+ |10101) - 100100} - |11001) - |00111) 
-100010} - 111100} - |00001) - |10000) 
-101110} - 110011} - 101000) + |11010) ] . 

The encoded Pauli operations for this code are conven- 
tionally expressed as 

X L = X ®X ®X ®X ®X 
Y L = Y ®Y ®Y ®Y ®Y 
Z L = Z®Z®Z®Z®Z. 

However, multiplying these encoded operations by mem- 
bers of the stabilizer group doesn't affect their action on 
the codespace. Thus we obtain the following equivalent 
set of encoded operations. 

X L = -X ® I ®Y ®Y ® I 
Y L = -Z ® Z <g> ï ®Y ® I 
Z L = -Y®Z®Y®I®I (13) 



These operators are all 3-local. This is the best that 
can be hoped for, because the code protects against 2- 
local operations and therefore any 2-local operation must 
anticommute with at least one of the generators. 

Besides increasing the locality of the encoded oper- 
ations, one can seek to decrease the number of qubits 
used to construct the codewords. The quantum single- 
ton bound|l7|| shows that the five qubit code is already 
optimal and cannot be improved in this respect. 

The distance d of a quantum code is the minimum 
number of qubits of a codeword which need to be modi- 
fied before obtaining a nonzero inner product with a dif- 
ferent codeword. For example, applying Xl, which is 
3-local, to |0l) of the 5-qubit code converts it into |1l), 
but applying any 2-local operator to any of the code- 
words yields something outside the codespace. Thus the 
distance of the 5-qubit code is 3. Similarly the distance 
of our 4-qubit code is 2. To detect t errors a code needs 
a distance of i + 1, and to correct t errors, it needs a 
distance of 2í + 1. 

The quantum singleton bound states that the distance 
of any quantum code which uses n qubits to encode k 
qubits will satisfy 

n-k>2(d-ï). (14) 

To detect 2 errors, a code must have distance 3. A code 
which encodes a single qubit with distance 3 must use at 
least 5 qubits, by equation^] Thus the 5-qubit code is 
optimal. To detect 1 error, a code must have distance 
2. A code which encodes a single qubit with distance 

2 must have at least 3 qubits, by eauation ll4l Thus it 
appears possible that our 4-qubit code is not optimal. 
However, no 3-qubit stabilizer code can detect all single- 
qubit errors, which we show as follows. 

The stabilizer group for a 3-qubit code would have two 
independent generators, each being a tensor product of 

3 Pauli operators. 

9i = cru ® 012 ® <Ti3 

g2 = CT 2 1 ® 022 ® a 23 

These must satisfy the following two conditions: (1) they 
commute, and (2) an X, Y, or Z on any of the three 
qubits anticommutes with at least one of the genera- 
tors. This is impossible, because condition (2) requircs 
o+: ^ o-ii 7^ I for each i = 1, 2, 3. In this case g\ and gi 
anticommute. 

The stabilizer formalism describes most but not all 
currently known quantum error correcting codes. We 
do not know whether a 3-qubit code which deteets all 
singlc-qubit errors while still maintaining 2-local encoded 
operations can be found by going outside the stabilizer 
formalism. It may also be interesting to investigate 
whether there exist computationally universal 3-local or 
2-local adiabatic Hamiltonians with a constant energy 
gap against local noise. 
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